Topological Entanglement Renyi Entropy and Reduced Density Matrix Structure 
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We generalize the topological entanglement entropy to a family of topological Renyi entropies 
parametrized by a parameter a, in an attempt to find new invariants for distinguishing topologically 
ordered phases. We show that, surprisingly, all topological Renyi entropies are the same, independent 
of a for all non-chiral topological phases. This independence shows that topologically ordered 
ground-state wavefunctions have reduced density matrices with a certain simple structure, and no 
additional universal information can be extracted from the entanglement spectrum. 



Introduction. Topological order (TO) [1] is a new 
kind of order that corresponds to patterns of long range 
quantum entanglement which cannot be described by 
symmetry breaking. However, the long range quantum 
entanglement in TO can leave its mark on the reduced 
density matrix, so one may be able to study long range 
entanglement and TO through the structure of these re- 
duced correlations. The reduced density matrix contains 
a lot of local non-universal information. The key is to fil- 
ter out all the non-universal information to capture the 
universal topological information, which is not affected 
by perturbations of the Hamiltonian, or small deforma- 
tions of the entanglement partition geometry. One way is 
to calculate topological entanglement entropy (EE) from 
reduced density matrices [2-6]. Such a universal quan- 
tity provides a way to determine whether or not a ground 
state possesses TO. If we only consider systems with a 
finite excitation gap, the low-energy physics can be de- 
scribed in terms of an underlying topological quantum 
field theory (TQFT). Then the topological EE is pro- 
portional to the logarithm of the total quantum dimen- 
sion Stop oc log-D. Unfortunately the quantum dimen- 
sion does not provide a complete classification of TO. For 
example, two topologically ordered states, the Z2 gauge 
theory and Ising anyons [7-9] , are different phases of mat- 
ter — with Abelian and non-Abelian anyonic excitations, 
respectively. However, they have the same Stop = log 2. 
To obtain a finer classification of TO, Ref. [10] proposes 
using the entire entanglement spectrum (possibly with 
additional conserved quantum numbers.) 

These developments motivated us to consider an ap- 
proach which might glean more universal information 
from the entanglement spectrum. We introduce a gen- 
eralization of the topological EE by deforming it into a 
Renyi entropy parameterized by a real number a which 
can characterize different aspects of the entanglement 
spectrum akin to moments of a probability distribution. 
We calculate this entropic quantity for the exactly solv- 
able string-net [11] and quantum double [12, 13] mod- 
els, which describe all the non-chiral topological phases. 
Recent works have mapped the quantum double models 



onto a subset of string-net models [14, 15], so we can 
compare entropies calculated for two different wavefunc- 
tions with the same TO. Our central result is that the 
only universal information captured by the Renyi entropy 
is the quantum dimension D, i.e. the topological Renyi 
entropy does not depend on the extra parameter a. As 
a consequence, no more universal information about the 
TO phases can be extracted from the entanglement spec- 
trum without additional conserved quantum numbers. 
Such a result suggests that the reduced density matrix 
Pa for a subregion A formally has the following structure 
®Pi = PA® Ptop, where ®pi is the tensor product of the 
local density matrices of the degrees of freedom living on 
the boundary of A. The "topological" desity matrix p tQV 
has a simple form where all its non-zero eigenvalues are 
equal, which leads to the a independence of the topolog- 
ical Renyi entropy, which we demonstrate explicitly for 
the quantum double models. 

Renyi entropy. The quantum Renyi entropy is de- 
fined with respect to a parameter a > as 



S a {p) 



1 



1 



log 



TV I 



(1) 



where the base of the logarithm is chosen to fix the units 
with which one measures the entropy. Taking the limit 
as a — > 1, one recovers the definition of the von Neu- 
mann entropy lim a _>i S a (p) = Si(p) = — Tr^plogp). 
The Renyi entropy is additive on independent states in 
the sense that the entropy of a product state is the sum 
of the individual entropies, S a {p <£> a) = S a {p) + S a (a). 
The Renyi entropy is essentially unique if we look for 
a function that is symmetric, continuous, has the addi- 
tive property, depends only on the spectrum of p, and 
obeys a generalized mean value property [16]. This es- 
sential uniqueness given certain natural assumptions and 
desired properties, together with the fact that the Renyi 
entropies cover a very broad class of functions motivates 
their consideration classification tool for TO. 

String-net states. We can study all parity-invariant 
topological phases in (2 + l)-d using string- net models 
[11]. These models exhibit TO and represent an exactly 
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solvable fixed point in a topological phase. The degrees of 
freedom are a set of strings living on the links of a hon- 
eycomb lattice. To specify a string net model requires 
several ingredients: a set of N string types i = 1, . . . , N, 
a branching rule tensor N % - k , and two real tensors di and 
^hi™ which satisfy certain algebraic relations [11] to en- 
sure consistency. Every string type i has an oppositely 
oriented partner i. The ground state wavefunctions of 
the string-net models obey a concise set of diagrammatic 
rules which are characterized by the string-net data listed 
above. In Ref. [4] the (von Neumann, a — > 1) topological 
EE for such string-net models was defined and calcu- 
lated to be <Stop = log D 2 where the quantum dimension 

m a knowledge of the ground 
state one can extract universal information about the 
low-energy TQFT and underlying TO in the form of the 
total quantum dimension. 

Renyi entropy for string-nets. To define an EE we 
begin by partitioning our system into two pieces. In this 
letter we will focus on a simply connected region A and 
trace out its exterior. The region A is topologically a 
disk and the reduced density operator of the string-net 
model on the disk can be deformed into a sum over string 
configurations on a tree-like diagram at the boundary of 
the disk [4]. We assume that our boundary string- net 
tree diagram has n boundary nodes with n links of the 
boundary tree labelled qi connected by n — 3 internal 
links. To begin the Renyi entropy calculation we start 
from Eq. 9 in Ref. [4], which gives the reduced density 
operator in region A, which we label by pa- We first raise 
PA to the power a and trace, summing over the states by 
using the branching rules iVj fe to get 
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where the expression for Ns q \ is given succinctly in 
terms of the matrices N q = £ a b \a)(b\, whose basis 
states form an orthonormal basis labelled by the string 
types: N {q} = (qi\N q2 N q3 ■ ■ ■ N qn _ 1 \q n ) . By relabelling 
the boundary strings in terms of the real- valued vector 
\d a ) = Y^,q dq \<l); we can return to Eq. (2) and write 
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where the sum on {q} runs only over n — 2 different qi. 
Since we are summing over all possible combinations, we 
can collect terms to get the even simpler form 
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We can make use of some properties of the N q matrices 
to simplify this expression. The N q satisfy iVj = N q 



(where q annihilates q) and if braiding is defined, we have 
N^ b = N£ , which implies that all the N q commute with 
each other. This means that the N q are normal and can 
be unitarily diagonalized simultaneously. Let S be the 
matrix such that S^NqS = A q is diagonal. Then we also 
have£ g JV? = S(E,A,dJ)St. 

Under the additional assumption that the braiding is 
sufficiently nontrivial (as discussed in the Appendix of 
Ref. [9]), we have so-called modularity, and the S de- 
scribed above is indeed the unitary modular S'-matrix of 
the theory. We choose the S-matrix to be in the canoni- 
cal form where we can read off the quantum dimensions 
from the first row or column. As we will see, this puts the 
largest eigenvalue of ~^2 q A q in the first matrix element. 

Since the N q are normal and mutually commuting, 
they share in common a complete set of orthogonal 
eigenvectors. Each N q has an eigenvalue d q with the 
eigenvector \d). Moreover, due to the Perron- Frobenius 
theorem, every other eigenvalue A for each N q satisfies 
|A| < d q . Thus we know exactly what the largest eigen- 
value of 'YliqNq is, namely ^ q d q . For symmetric ma- 
trices (and y_2 q N q is symmetric), the Perron- Frobenius 
theorem gives us additional guarantees. In particular, 
the largest eigenvalue A max is non-degenerate. Further- 
more, the least eigenvalue satisfies A m i n = — A max if and 
only if the symmetric matrix is the adjacency matrix of 
a bipartite graph. But this can't be the case, since the 
vacuum always fuses with itself to form the vacuum, giv- 
ing at least one nonzero element on the main diagonal, 
and bipartite graphs have no self-loops. Therefore all 
other eigenvalues A of J2 q N q satisfy |A| < J2 q d q , and 
these A contribute exponentially less once we raise to the 
power Ti — 2. Then, ignoring a multiplicative factor of 
(1 + 0(cxp(— n))), we have 
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To get a more explicit expression, we need to calculate 
{cZ a 1) | 2 . Let's consider how S acts on |1). S is a 
unitary matrix, and the first row is proportional to (d\. 
So S\l) = jg\d). Hence |(d Q |S*|l)| 2 = (d a \d) 2 /D 1 and 
substituting this in to Eq. (5) and using the expression 
for the Renyi entropy in Eq. (1), we obtain 
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which is correct up to a term of order 0(exp(— n)). The 
first term represents the area law. It is not universal and 
cannot be used to describe the phases. The second term 
represents the universal part: the topological entangle- 
ment Renyi entropy. We see that it does not contain 
any a dependence, just the total quantum dimension D. 
Therefore the it does not provide any additional universal 
information. The Renyi entropies completely determine 
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the spectrum, hence no additional information (beyond 
D) can be gathered from the entanglement spectrum. 
This is true when the partition geometry is simply con- 
nected; Ref. [17] has shown that more can be extracted 
in more complicated partitions. 

We wish to find deeper insight into why there is noth- 
ing else in the eigenvalues of the reduced density ma- 
trix that can say more about topological order. To this 
end, we will study an important class of TO states, those 
emerging from discrete gauge theories. In the following, 
we prove that the reduced density matrix of such states is 
proportional to a projector, and thus all Renyi entropies 
contain no a dependence and that the whole entangle- 
ment spectrum is trivial and flat. 

Quantum Double Models. The quantum double mod- 
els are exactly solvable lattice models with discrete gauge 
symmetries [12, 13]. These models exhibit phases with 
TO and anyonic excitations, and are in the same univer- 
sality class as a subset of the string-net models [14]. To 
define them, begin with a directed graph with orienta- 
tions ± and with qudits on the edges. Consider a finite 
group G of dimension |G| = d, with identity e. The local 
Hilbcrt space on the edge i is therefore Hi ~ C[G] and an 
orthonormal basis for the qudits is given by {\g) : g G G}. 
The total Hilbert space for a system with n qudits is given 
by H = <3™ =1 Hi. We focus on the model on a square lat- 
tice, with n/2 vertices and plaquettes. 

Following the construction of [13], the relevant op- 
erators are L 9 ±,T£ defined by L 9 + \z) = \gz),T^\z) = 
5 htM \z),LL\z) = \zg~\Th\z) = 5 h -i,,\z). The gauge 
transformations are defined as follows: 

4 

M*) = I[L 9 (j,s) , MP) = E II Thm U^P)- 

j£s hxh^h^h^— e m— 1 

The star and plaquette operators are defined as the 
projector operators A(s) = | C | 1 J2 g( =G ^s( s )' B(p) = 
B e (p). The Hamiltonian of the quantum double model is 

h qd = y,( 1 - a (s)) + Y,( 1 - b (p)) ( ? ) 

s p 

Since [A(s),B( P )} = [A(s),A(s')] = [B(p),B(p')} for all 
s,s',p,p', the ground state manifold is given by the set 
£ = {|0 € H | A(s)\0 = B(p)\0 = |£) Vs,p} with 
ground state energy Eq = 0. 

Consider the vacuum state |e) = |e)®". For each pla- 
quette p, it easily follows that B(p)\e) = \e). We can 
build a (un- normalized) ground state |£o) € C by pro- 
jecting as follows, |£o) = Els A(s)\e). 

Now, consider the set ® of all the possible A g {s). What 
this operator does is to make a small loop around s with 
string of type g. We have = {A g (s),g G G, s = 
1, ...,n/2}. Now consider the set Q = (©}, that is the 
set of all the possible products of elements in 25. The set 



Q is a group. With this definition, we have 

i&> = \g\- x n E A ^)V) = i G r f E %> ( g ) 

s geG fi6S 

It is important to see that the set of {\h)} is orthonormal. 
Moreover, given a bipartition of the Hilbert space % = 
Ha ®H-b, the set <8> \hs)} is bi-orthonormal. Let 

us compute the density matrix po = |£o)(£o|- Since each 
vector \h) factorizes as \Ka) 8> |/ib), we have 

Po = \G\- n E \h A ){h' A \®\h B )(h' B \, (9) 
hM'eg 

Consider now the subgroup of Q acting exclusively on 
subsystem A, Qa := {g G Q \ g = gA <8 is}, 
and analogously consider Qb- It is easy to show that 
Qa, Qb, and Qa x Qb are normal in Q. Therefore 
we can define the quotient groups Qab '■= Q/Qa x 
Qb,Q/Qb,Q/Qa- We see that the only elements of Q 
such that (e\fiB\e) ^ are those in Qa, and therefore we 
find p A = \G\~ n J2h&g,h&SA I^a)^ 1 ^!) wner e we have 
relabelled the group elements as h' = hr x h. Notice that 
\h A )(hA\ = g\h A )(h A \ for every g G Qb, and \Q\ = \G\ n . 
Therefore, reordering gives 

pa = \Q\~ 1 \Qb\ E lO(^M (io) 

heg/9B,heQA 

Squaring this expression for pa and using the group 
properties shows that pa is proportional to a projector, 
Pa = ^ A \g\ B ^ PA, and therefore the Renyi entropies con- 
tain no additional information beyond the quantum di- 
mension, D = \G\. The entanglement spectrum is flat 
which is connected [10] with the trivial nature of the edge 
states for the QD models. 

The origin of the topological term. At this point, we 
would like to understand why pa is just a projector? And 
why, in the more general string-net setting where the re- 
duced density matrix is not just a projector, is there still 
no topological information other than D? Here we prove 
that the reduced density matrix pa is unitarily equiva- 
lent to a matrix that only addresses the degrees of free- 
dom on the boundary of the partition. Moreover, we 
show that the area law has a correction because there 
is a global constraint on the boundary. We can enlarge 
the system by removing this constraint and express the 
reduced density matrix as the tensor product of the lo- 
cal density matrix of each of the degrees of freedom on 
the boundary We focus on the Z2 case for simplicity, 
but the argument can be generalized to all the quantum 
double models. In this case, the ground state is given by 
Eq. (8), where Q is the group generated by the plaque- 
tte operators A v = Iljedp^f and |0) is the state with 
all spins up in the z-basis. By choosing a simply con- 
nected region of plaquettes, we partition the spins into 
{A, £>), where A includes the spins in the interior and on 



4 



the boundary. The quotient group Qab consists of the 
closed strings that act on both A and B, that are equiv- 
alent under deformations acting entirely within A or B. 
Therefore, the equivalence classes in Qab can be repre- 
sented by those closed strings that live near the bound- 
ary between A and B, namely those closed strings that 
are generated by the plaquettes that are external to A 
and share one edge with the boundary. So every element 
h G Gab can be decomposed as h = Ka ® h B where Ha 
only acts on spins that live on the boundary ( "a" spins) . 
The h B part only acts on those spins which are external 
to A ( "6" spins) . The rest of the lattice consists of the 
spins in the bulk of A and B, namely all those spins that 
belong solely to either A or B: |0) = |0) a <g> |0) 6 <g> |0) bu ik 
so that h\0) = h A \0) a <g> h B \0) b ® /|0) bu ik- Therefore the 
ground state can be written as 



gA®gB£0AxQB 
heQAB 



h A 1 0) a ® h B 1 0) b <g> (g A <8> 9b ) | 0) 



bulk 



\Q 



X 



-1/2 
1-1/2 



gx, with X = A,B. We 



\gAB\- 1/2 QAQBEh & g AB h A\0)a®h B \0}b® 



Define Qx 
obtain \ip) - 

|0)b u ik- The density matrix can be therefore be factored 
as p = QaQbP® P (hulk) QAQB, with 

p®p< bulk ) = \Qab\~ 1 h A\Q)(V\ati A ®h B \0)(0\bti B 



®|0)<0|buik. 



Notice that the bulk part is separable in the bipartition 
(A, B) so that the reduced density matrix can be written 

as pa = QaPa <8> Pa^Qa, with p^ nlk ^ being a pure 
state, so the EE of p is just the entropy of pa- Then 



PA = \Q, 



AB 



^2 h A \0)(0\ a h A) 



(11) 



which gives the expected result for the entanglement 
Si(p) = log 2 |^AB| [2, 18]. We are interested in un- 
derstanding what are the spin configurations in the sum 
Eq. (11). Notice that the support of pa consists only of 
the spins on the boundary. Now, note that every spin 
configuration is not allowed. In fact, we have the follow- 
ing global constraint: Y\ he g AB h = g A ®> gs G Qa x Qb- 
So the product of all the Ha is also in Qa and is +1 
on the ground state. The global constraint is thus 
YijedA a j = +1 so that the reduced density matrix 
Eq. (11) consists of the sum of all spin configurations 
with parity +1, namely all the spin configurations with 
an even pair of spins flipped. If the boundary has length 
n, then there are 2™ _1 such configurations and the entan- 
glement is then S = n — 1. Now we understand that the 
topological state is completely determined by the bound- 
ary, and that we have the completely mixed state within 
the sector of parity +1. We can consider the enlarged 



system by considering the perfect mixture with the sec- 
tor of parity — 1. In this case, we have that 



(area) 

Pa 



PA 



1 
1 



- _ (top) 

PA® Pa ■ 



(12) 



So Pa is just the completely mixed state of all the 
possible spin configurations on the boundary and thus 



(area) 

Pa 



®i=iPi 



- _ (top) 

PA® P A ■ 



(13) 



We have shown that the entanglement in the ground state 
of a topologically ordered system is completely contained 
in the boundary, namely in the entropy of the reduced 
density matrix p~A- We have also shown that this state 
almost obeys an area law, because there is a global topo- 
logical constraint, namely that only spin configurations 
of parity +1 are allowed. Therefore, we can complete 
it with a density matrix that describes a system before 
we project onto the system with parity +1 . This term 
contains the topological entropy. Once completed, the 
system obeys a strict area law and decomposes into the 
local tensor product of the single degrees of freedom on 
the boundary. Such structure of the reduced density ma- 
trix, as described in Eq. (13), explains why the topolog- 
ical Renyi entropies do not depend on a. 

Finally, we remark that our proof applies to non-chiral 
topological phases. Therefore it is still an open problem 
to what extent the entanglement spectrum can classify 
chiral topological phases [10, 19]. 
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